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ABSTRACT 

In the course of its commissioning observations, the Sloan Digital Sky Survey (SDSS) has 
produced one of the largest redshift samples of galaxies selected from CCD images. Using 11,275 
galaxies complete to r* = 17.6 over 140 square degrees, wc compute the luminosity function of 
galaxies in the r* band over a range —23 < A'/,-* < —16 (for h = 1). The result is well-described 
by a Schechter function with parameters (f)^ = (1.46±0.12) x 10"^ Mpc'^, = -20.83±0.03, 
and a = -1.20 ± 0.03. The implied luminosity density in r* is j Ri (2.6 ± 0.3) x lO^/iLoMpc"^. 
We find that the surface brightness selection threshold has a negligible impact for M,.* < —18. 
Using subsets of the data, we measure the luminosity function in the u*, g*, i*, and z* bands 
as well; the slope at low luminosities ranges from a = —1.35 to a = —1.2. We measure the 
bivariate distribution of r* luminosity with half-light surface brightness, intrinsic g* — r* color, 
and morphology. In agreement with previous studies, wc find that high surface brightness, red, 
highly concentrated galaxies are on average more luminous than low surface brightness, blue, 
less concentrated galaxies. An important feature of the SDSS luminosity function is the use of 
Petrosian magnitudes, which measure a constant fraction of a galaxy's total light regardless of 
the amplitude of its surface brightness profile. If we synthesize results for i?GKC-band or fe^-band 
using these Petrosian magnitudes, we obtain luminosity densities 2 times that found by the Las 
Campanas Redshift Survey in i?GKC and 1.4 times that found by the Two-degree Field Galaxy 
Redshift Survey in bj. However, we are able to reproduce the luminosity functions obtained by 
these surveys if we also mimic their isophotal limits for defining galaxy magnitudes, which are 
shallower and more redshift dependent than the Petrosian magnitudes used by the SDSS. 

Subject headings: galaxies: fundamental parameters — galaxies: photometry — galaxies: statis- 
tics 



1. Motivation 

A fundamental characteristic of the galaxy population, which has been the subject of study at least 
since Hubble (1936), is the distribution of their luminosities. Given the broad spectral energy distributions 
of galaxies (especially in the presence of dust), the bolometric luminosity of galaxies is currently too difficult 
to measure to meaningfully study, so we must be satisfied with observing the luminosity function in some 
wavelength bandpass. In this paper wc measure the luminosity function of local {z < 0.2) galaxies in five 
optical bandpasses, between about 3000 A and 10000 A. We also investigate the correlation of luminosity 
with other galaxy properties, such as surface brightness, intrinsic color, and morphology. This quantitative 
characterization of the local galaxy population provides the basic data that a theory of galaxy formation 
must account for and an essential baseline for studies of galaxy evolution at higher redshifts. 

The most recent determinations of the optical luminosity function of "field" galaxies (those selected 
without regard to local density) have come from large flux-limited redshift surveys. They include the lu- 
minosity function measurements of Lin et al. (1996), using the Las Campanas Redshift Survey (LCRS; 
Shectman et al. 1996) in the LCRS iicKC-band (around 6500 A), and of Folkes et al. (1999) using the 
Two-degree Field Galaxy Redshift Survey (2dFGRS; CoUess 1999) in the b^-band (aroimd 4500 A). Other 
recent determinations of the local optical luminosity function have been mostly in the B- or 6^ -bands, for 
example the Nearby Optical Galaxy sample (Marinoni et al. 1999), the Optical Redshift Survey (Santiago 
et al. 1996), the Stromlo-APM Redshift Survey (Loveday et al. 1992), the Durham/UKST Galaxy Redshift 
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survey (RatclifFe et al. 1998), and the ESO Slice Project Galaxy Redshift Survey (Zucca et al. 1997). Of 
particular interest here is the ability to use the multi-band photometry of the SDSS to meaningfully compare 
our results to those of these other surveys. 

In this paper, we use imaging and spectroscopy of 11,275 galaxies over about 140 square degrees complete 
to r* = 17.6 from the Sloan Digital Sky Survey (SDSS; York et al. 2000) commissioning data to calculate 
the luminosity function of galaxies. In Section 2 we describe the source catalog and the redshift sample, 
with particular attention to our definition of galaxy magnitudes, which is based on a modified form of the 
Petrosian (1976) system. In Section 3, we describe the maximum-likelihood method used to calculate the 
luminosity function. In Section 4, we present the results for the luminosity function in each band, for the 
luminosity density of the universe, and for the dependence of galaxy luminosity on surface brightness, color 
and morphology. In Section 5, we examine the effects of using Petrosian magnitudes rather than isophotal 
magnitudes for the luminosity function. In Section 6, we compare our results to those of Lin et al. (1996; 
LCRS) and Folkes et al. (1999; 2dFGRS). We show that the SDSS luminosity function implies a substantially 
higher mean luminosity density than cither of these previous measiircmcnts, hut that we can reproduce the 
results of both surveys if we mimic their isophotal definitions of galaxy magnitudes. We conclude in Section 
7. 



2. SDSS Commissioning Data 

2.1. Description of the Survey 

The SDSS (York et al. 2000) will produce imaging and spectroscopic surveys over tt steradians in the 
Northern Galactic Cap. A dedicated 2.5m telescope (Siegmund et al., in preparation) at Apache Point 
Observatory, Sunspot, New Mexico, will image the sky in five bands (u', g', r', i', z'; centered at 3540 A, 
4770 A, 6230 A, 7630 A, and 9130 A respectively; Fukugita et al. 1996) using a drift-scanning, mosaic CCD 
camera (Gunn et al. 1998), detecting objects to a flux limit of r' 23. Approximately 900,000 galaxies, 
(down to rfjjjj w 17.65), 100,000 Bright Red Galaxies (BRGs; Eiscnstein et ai, in preparation), and 100,000 
QSOs (Fan 1999; Newberg et al., in preparation) will be targeted for spectroscopic follow up using two 
digital spectrographs on the same telescope. The survey has completed its commissioning phase, during 
which the data analyzed here were obtained. 

As of January 2001, the SDSS has imaged around 2,000 square degrees of sky and taken spectra of 
approximately 100,000 objects. In this paper, we concentrate only on ^ 230 square degrees along the 
Celestial Equator in the region bounded by 145° < a < 236° and -1.25° <S < 1.25° (J2000). This region 
was imaged during two runs, known as SDSS runs 752 and 756, over two nights in March 1999. Each run 
consists of six columns of data, each slightly wider than 0.2 degrees and separated by about the same amount; 
the nms are interleaved to form a complete, 2.5 degree wide, stripe. The seeing varied over the course of 
the runs from about 1.2" to 2", with a median of approximately 1.5". In this area, the spectroscopic survey 
has so far obtained more than 11,275 redshifts of galaxies selected from these runs, over about 140 square 
degrees of sky, from which we here calculate the luminosity function. 
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2.2. Identifying and Measuring Galaxies 

The images produced by the camera are analyzed by a photometric pipeHne specifically written for 
reducing SDSS data (photo; Lupton et al, in preparation), photo subtracts the CCD bias, divides by the 
flat field, interpolates over data defects (cosmic rays, bad columns, bleed trails), finds objects, dcblends 
overlapping objects, and performs photometry. It estimates the local point-spread function (PSF) as a 
function of position based on bright stars. Among many other tasks, for each photometric band the pipeline 
calculates the 3" diameter aperture magnitude mgber, the PSF magnitude mpsF (that is, the magnitude 
using the local PSF as a weighted aperture), and the Petrosian magnitude mp (described more fully below). 
It fits for the parameters of (possibly elliptical) de Vaucouleurs and exponential profiles for each object. It 
picks the better model of the two in r* , and uses this model as a weighted aperture to calculate the "model 
magnitude" mmodci in sXl bands; as described in the Section 2.4, while the fiux limit is defined with repect 
to Petrosian magnitudes, the model magnitudes are used to distinguish between stars and galaxies. More 
sophicated models will be fit to a bright subsample of the galaxies in a later paper, but the survey pipelines 
measure model magnitudes for every detected object in the survey, and to fit more complicated models to 
all objects would be prohibitively expensive. For each object, photo also measures an azimuthally averaged 
radial profile out to the maximum detectable distance. 

The measured magnitudes are calibrated in the ABg^ system (Fukugita et al. 1996). The system is 
based on three "fundamental standards" (BD +17°4708, BD +26°2606, and BD +21°609). The United 
States Naval Observatory Im telescope has calibrated about 150 "primary standards" with respect to the 
fundamental standards. Meanwhile, a 0.5m photometric telescope (PT; Uomoto et al., in preparation) is 
used to relate the primary standards to a set of secondary standards, all fainter than 14™ and lying in patches 
referred to as "secondary calibration patches." These patches lie along the survey stripes. Thus, during the 
course of its imaging on any given night, the c;anicra passes ovcir these secondary calibration patches; because 
the secondary standards are faint enough that they do not saturate the camera on the 2.5m, they can be 
used to calibrate yet fainter objects detected by the imaging. Then, all of the detected objects can be related 
to the primary standards, and thus to the fundamental standards. The PT serves yet another function, 
which is to monitor the primary standards during nights that the 2.5m telescope is imaging, in order to 
determine extinction coefficients and their variation over the course of the night. While this procedure was 
used to calibrate the commissioning data, the resulting calibration has not yet been fully verified. Here, we 
represent magnitudes with the notation u*, g*, r* , i*, and z*, rather than the standard notation u', g', r', 
i' , and z\ to indicate the preliminary nature of the calibrations used here. We expect the final system to 
differ from the current system by 0.05 magnitudes or less. In particular, we have found consistent results 
from calculating the luminosity function in regions other than the equatorial stripe presented here, which 
were calibrated based on different secondary patches. 



2.3. Petrosian Magnitudes 

Because galaxies are resolved objects with poorly defined edges, and do not all have the same radial 
surface brightness profile, some care is required to define the "flux" associated with each object. An ideal 
method would measure the "total" light associated with each galaxy, but in practice any such method 
requires a model-dependent extrapolation of the measured light profile, and an accurate extrapolation itself 
requires high signal-to-noise ratio data in the outer regions of the galaxy. A conventional alternative is to 
measure flux within a specified isophote, but with such isophotal magnitudes the fraction of a galaxy's light 
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that is measured depends on the amplitude of its surface brightness profile, and the fraction decreases if the 
surface brightness is diminished by cosmological redshift dimming or by Galactic extinction. To avoid these 
problems, the SDSS has adopted a modified form of the Petrosian (1976) system, measuring galaxy fluxes 
with a circular aperture whose radius is defined by the shape of the galaxy light profile. 

More specifically, we define the "Petrosian ratio" TZp at a radius r from the center of an object to be 
the ratio of the local surface brightness averaged over an annulus at r to the mean surface brightness within 
r: 



where /(r) is the azimuthally averaged surface brightness profile and aio < 1, ahi > 1 define the annulus. 
The SDSS has adopted aio = 0.8 and am = 1.25. 

The Petrosian radius rp is defined as the radius at which TZp{rp) equals some specified value T^p^um- 
The Petrosian fiux in any band is then defined as the flux within a certain number Np of r* Petrosian radii: 



Thus, the aperture in all bands is set by the profile of the galaxy in r* alone. The SDSS has selected 
T^p.iim = 0.2 and Np = 2. The aperture 2rp is large enough to contain nearly all of the light for a typical 
galaxy profile (see below), so even substantial errors in rp cause only small errors in the Petrosian flux, but 
small enough that sky noise in Fp is small (typical statistical errors near the fiux limit of r* = 17.65 are 
< 5%). In practice, there are a number of complications associated with this definition: the galaxy images 
are pixelized, some galaxies have undefined Petrosian radii because sky noise begins to dominate before 
the Petrosian ratio drops to T^p, some galaxies have multiply-defined Petrosian radii because of galaxy 
substructure, and so forth. We defer detailed discussion of these issues to Lupton et al. (in preparation), 
since for nearly all galaxies in the spectroscopic sample the idealized account above is accurate. 

Given the Petrosian flux, one can find the radius rso contained 50% of the Petrosian flux and the radius 
rgo containing 90% of the Petrosian flux. As shown below, the Petrosian flux measures virtually all of the flux 
in pure exponential profiles but only 80% of the fiux in pure de Vaucouleurs profiles. Thus, for exponential 
profiles, r5o and rgo correspond to the true 50% and 90% light radii, but for de Vaucouleurs profiles r5o,petro ~ 
0.7r5o.truc and rgo,pctro ~ 0.43rgo,truc- We define the half-light surface brightness to be the average surface 
brightness within the half-light radius, in magnitudes per square arcsecond: /Xi/2 = inp + 2.5 log]^Q(27rr|Q). 
The "concentration index" of galaxies is defined as c = rgo/r5o. High concentration index objects, such 
as de Vaucouleurs profile galaxies, have a strong central concentration of light and large, faint wings. Low 
concentration index objects, such as exponential profile galaxies, have a light distribution which is closer to 
uniform. We will generally express our results as a function of 1/c, the inverse concentration index. 

The Petrosian ratio Tip \s manifestly indifferent to multiplicative changes in the surface brightness of 
the galaxy and depends only on the physical radius in the galaxy, not on its redshift or Galactic extinction. 
Thus, the ratio of the Petrosian flux to total flux for any given galaxy depends on redshift only to the extent 
that the effects of seeing and i^-corrections change the measured shape of the galaxy profile. Petrosian 
magnitudes only become practical for relatively deep imaging; otherwise the Petrosian ratio (Equation 1) is 
noisy and the aperture cannot be made large enough to capture most of a typical galaxy's light. 

Using simulated galaxy observations, we can determine what fraction of the light the Petrosian flux 
contains for typical galaxy profiles. This fraction is independent of redshift except when the size of the 




(1) 




(2) 
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galaxy is comparable to the seeing. Figure 1 shows the ratio Fp /Ftotai of the Petrosian flux to the total flux, 
as a function of the (true) angular half-light radius of a galaxy. The top panel shows Fp /i^totai for a face- 
on exponential disk; the bottom panel shows the same quantity for a circularly symmetric de Vaucouleurs 
profile. The dotted line is the result in the limit of no seeing, for an axisymmetric galaxy. The solid line 
is the result assuming the median seeing of 1.5" for the runs analyzed here, again for an axisymmetric 
galaxy. To calculate these quantities, we convolve each model galaxy with a realistic point-spread function 
(including the power-law wings) measured from bright stars in the survey. As the size of the galaxy becomes 
comparable to the seeing, the fraction of light enclosed within a radius of 2rp becomes closer and closer to 
the corresponding fraction for a star, which is about 95%. In the case of exponential profiles, Fp « -Ftotai m 
the absence of seeing. Seeing therefore reduces the fraction of light measured by the Petrosian flux, though 
only by about 3% even for small galaxies. In the case of de Vaucouleurs profiles, about 82% of the light 
is measured by the Petrosian flux in the absence of seeing. Seeing therefore increases the fraction of light 
measured for a small de Vaucouleurs galaxy by about 10%. These results depend very little on the inclination 
of the disk or the axis ratio of the de Vaucouleurs profile. For example, in Figure 1, the dashed line in each 
panel shows the result for a galaxy with an axis ratio of 0.5, convolved with seeing. 

For comparison, we show in Figure 2 the distribution of observed half-light radii. This figure shows the 
distribution of r^Q for galaxies with r* < 17.6 in the top panel, and the fractional cumulative distribution 
-^(< fdo) in the bottom panel. Roughly half of the galaxies have r^o < 2.5", and thus probably have their 
Petrosian fluxes under- or overestimated by a few percent. We emphasize that a comparably deep isophotal 
magnitude would have a similar dependence on angular size to that found here, without the benefits of being 
independent of galaxy redshift for galaxies with half-light radii larger than the seeing. 

2.4. Targeting Galcixies 

Once objects have been identified and photometry has been performed, which can be done reliably to 
r* ~ 23, a sample of targets is selected for spectroscopic observation. The details of this selection and the 
software pipeline used to implement it will be described in separate papers (Strauss et al., in preparation; 
vanden Berk et ai, in preparation). For the purposes of the main galaxy sample, extended objects (galaxies) 
are separated from stars based on the quantity mpsp — Wmodei in where (as described in Section 2.2) 
mpsF is the magnitude using a PSF-weighted aperture and mmodei is the magnitude associated with the 
best-fit model to the galaxy profile. That is, if most of the object's light is contained within an aperture 
weighted by the PSF, it is considered to be stellar. Figure 3 shows the distribution of our spectroscopic 
galaxy targets versus mpsF — "imodei- The objects which are galaxies (based on the spectroscopic data) 
are shown in the solid histogram, and those which are stars are shown in the dotted histogram. For the 
commissioning data, the division was set at 0.242; roughly 1% of objects brighter than the spectroscopic 
limit and identified as galaxies in this way are stars. Indeed, most of the stellar spectra recovered from the 
galaxy sample are actually double stars, which shows that this method of separating the populations is fairly 
clean. It is clear from the drop-off in the distribution of galaxies as mpsF — Wmodei decreases that galaxies in 
this magnitude range with mpsp — Wmodci < 0.242 are extremely rare. We use the model magnitude here for 
star/galaxy separation rather than the Petrosian magnitude because the model magnitude is more robust at 
faint magnitudes (r* > 21) and to the effects of seeing; we emphasize that for the spectroscopic sample such 
considerations are irrelevant, because the star /galaxy separation is so clean. 

Given the set of extended objects in the commissioning data, wc chose as spectroscopic targets all objects 
with a Petrosian magnitude (after an extinction correction using the dust maps of Schlegel, Finkbeiner & 
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Davis 1998) of r* < 17.65 and a fraction of those between 17.65 and 17.75, with a probability dechning 
linearly from unity to zero in this range. This flux limit was applied using the calibration available in the 
March 2000, in a version of the processing known as "rerun 1." The magnitudes quoted here were calibrated 
in a later version known as "rerun 4." Because of the differences in calibration, to define a complete sample 
we have restricted ourselves here to a flux limit of r* < 17.6. 

In addition, a surface brightness cut of /ix/2 < 23.5 was used for the commissioning observations; this 
cut eliminates less than 1% of all extended objects brighter than the flux limit, about half of which are 
real galaxies and half of which are ghosts or consequences of improper deblending. We have checked the 
effects of improper deblending by examining by eye all of the targets in the range 22.5 < /Ui/2 < 23.5. 
Fifteen (about 3%) of the targets in this range arc products of deblending failures which nevertheless have 
galactic redshifts successfully determined and thus contaminate our sample; we have checked that the effect 
of this contamination is small compared to our estimated errors. The main survey uses improved methods 
for eliminating contamination by such defects and will accordingly have a somewhat more relaxed surface 
brightness limit. We will explore below (Section 4.3) some of the consequences of the surface brightness limit 
in the commissioning data. 

A number of other requirements were used in selecting spectroscopic targets. The "fiber magnitude" 
(defined in Section 2.2 as the magnitude within an 3" diameter aperture) was restricted to be fainter than 
15 in g*, r*, or i* to avoid cross-talk between fibers on the spectrograph CCDs. Objects brighter than 
r* = 15.5 and with half-light radii < 2" almost invariably turn out to be misclassificd stars and therefore 
we reject them; this cut always affects < 1% of objects, and only approaches this percentage in poor seeing 
conditions. We note that the BRG sample targets galaxies fainter than r* = 17.65 that satisfy specified 
color and magnitude selection criteria; we have not made use of the BRG sample in the luminosity function 
measurements presented here. 

2.5. Spectroscopic Observations 

The spectra of the chosen galaxies are obtained using a multi-object fiber spectrograph that can take 
640 spectra at once in a circular plate (or "tile") 1.49 degrees in radius (Uomoto et al, in preparation). 
Each exposure is about 45 — 60 minutes, depending on observing conditions. Each fiber has a diameter of 
3", and is led into a dichroic beam-splitter centered at ~ 6000 A which separates the light onto red and blue 
cameras. The resulting spectrum covers the range from 3800 A to 9200 A, with a resolution R 2000. 

In practice, 48 of the fibers are used as sky fibers and spectrophotometric standards, leaving 592 for 
quasars, galaxies, and stars. Typically, about 500 fibers per plate are used for galaxies. Although this 

number provides enough fibers for all our desired targets, the galaxies arc distributed non-uniformly on the 
sky. Furthermore, no two fibers on the same plate can be placed closer than 55", a limit imposed by the 
physical size of the fiber plugs. For these two reasons, we must position the centers of the spectroscopic 
plates and allocate the targets to the plates carefully, in order to maximize the number of targeted galaxies 
(Blanton et at, in preparation). In the commissioning data we assigned fibers to about 92% of the targeted 
objects (lumping main sample galaxies, QSOs, and BRGs together), although for the survey proper, for 
which the geometry will be considerably more favorable, our sampling will be above 95%. Most of the 
missing targets are due to the 55" collision constraints. 

After observation, the spectra are analyzed using a specially written spectroscopic pipeline (Prieman 
et al., in preparation) which extracts the 640 spectra from each two-dimensional spectrogram, carries out 
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wavelength and flux calibration, subtracts the sky, classifies the objects as stars, QSOs, or galaxies, and 
determines the redshift. For galaxies, redshifts are determined independently from emission lines and from 
a cross-correlation with absorption spectra; using both methods provides a useful check on the software. In 
the data analyzed here (not all of which is considered "survey quality"), and with a commissioning version of 
the pipeline, the redshift success rate for objects targeted as galaxies is about 98%, as determined by eyeball 
examination of the spectra, independent of apparent magnitude and surface brightness. The remaining 2% 
fail to give redshifts, due for the most part to broken fibers and the accidental targetting of ghost images. 
The redshift errors are around 30 km s~^, as determined by repeat observations of the same objects; this 
error is completely negligible for the purposes of measuring the luminosity function. 



2.6. Sampling Rate 

To calculate the luminosity function, we must map the selection function of the survey. Aside from the 
flux limit, there are three main effects: missing galaxies due to lack of flbers in dense regions, missing galaxies 

due to spectroscopic failures, and missing galaxies due to fiber collisions. To account for these effects, we 
calculate the local sampling rate of galaxies separately for each region covered by a unique set of plates. 
That is, in the case of two overlapping plates, the sampling is calculated separately in three regions: the 
region covered only by the first plate, the region covered only by the second plate, and the region covered 
by both. We refer to each of these three regions as an "overlap region." Thus, depending on what overlap 
region it was in, each galaxy can be assigned a sampling rate ft, equal to the number of redshifts of galaxy 
targets obtained in the region divided by the number of galaxy targets in the region. In regions covered by a 
single plate, ft ~ 0.85-0.9; in multiple plate regions, ft ~ 0.95-1. These completenesses average to the 92% 
quoted above. 

It is possible to simply weight each galaxy by 1 /ft when calculating the luminosity function. However, 
the majority of the missing galaxies are missing because of the fiber collision constraints. It is more appropri- 
ate to account for the collisions as follows. First, we group the galaxies using a friends-of-friends algorithm 
with a 55" linking length. For each galaxy for which we obtained a redshift, we assign a "collision weight" 
Wc equal to the multiplicity of the group divided by the mimbcr of redshifts obtained in the group. That is, 
in a colliding pair for which we have one redshift, we weight that redshift by a factor Wc = 2. This procedure 
is essentially the same as assuming the two galaxies have the same absolute magnitude. The sampling rate 
using these collision weights, which wc denote ft, turns out to be ^ 0.97 1 for most overlap regions. When 
we calculate the normalization of the luminosity function below, we weight each galaxy by the factor wd ft- 
We should note that the difference between weighting by wd ft and by 1/ft is negligible in the calculation 
of the luminosity function (although this is certainly not the case for clustering statistics). 



2.7. Distance Modulus and ii'-corrections 

To convert the apparent magnitude m to an absolute magnitude M, we must assume a particular 
cosmology and correct for the fact that the observed bandpass differs from the restframe bandpass. The 
relation between m and M can be written 

M = m- DM(z) - K{z), (3) 



where DM(z) is the bolometric distance modulus and K{z) is the iiT-correction. 
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Here we use distance moduli based on three different Friedman- Robertson- Walker models: {ilm = 1, 
f^A = 0); {ilm = 0.3, = 0.7); and (fJ^ = 0.3, = 0). See Hogg (1999) and references therein 
for distance measure formulae for these cosmologies. At z = 0.2, the highest redshift considered here, the 
difference in the distance modulus for these three cases is about 0.2 magnitudes at most. Meanwhile, the 
comoving volume associated with a coordinate volume clement dfldz differs among these models by as much 
as 40%. We use all three models for the calculation of our main results for the luminosity in each band, but 
for the investigation of the dependence of luminosity on other parameters, we consider only the Clm = 0.3, 
f^A = 0.7 model for the sake of simplicity. We use Hq = 100ft. km s~^ Mpc~^, with h = 1 when we have 
not specified h explicitly; choosing alternative values for h shifts all absolute magnitudes by 5 log^g ^^'^ 
changes number densities by the factor h^. 

Because we observe the galaxies in a fixed observed wavelength band, a /sT-correction is necessary to 
account for the fact that the observed band corresponds to different rest-frame bands at different redshifts. 
Naturally, this ii'-correction is dependent on the galaxy spectral energy distribution (SED). Given the five- 
band SDSS photometry and the high quality spectra over a large range in wavelength, the best approach in 
the long run will be to determine the iiT-corrections of the galaxies from the survey itself. For the present, 
however, we will simply use ii'-corrections supplied by Fukugita, Shimasaku, & Ichikawa (1995) for a range 
of galaxy types. Using the observed galaxy redshift and the g* — r* color, we interpolate between the galaxy 
types in Fukugita, Shimasaku, & Ichikawa (1995) to determine the intrinsic g* — r* color and the appropriate 
iC-correction. Figure 4 shows the iiT-correction in r* as a function of redshift for several different galaxy 
colors over the range of redshifts of interest here. While the true survey bandpasses are not exactly those 
assumed by Fukugita, Shimasaku, & Ichikawa (1995), we have tested the effect of using ii'-corrections based 
on recent estimates of the filter curves (Fan et al. 2001), and the differences in the resulting measured 
luminosity density are less than 5% in all bands. 



3. Calculating the Luminosity Function 

We define the unit-normalized luminosity function $(L) to be the distribution function of galaxies in 
luminosity, normalized to unit integral over the range from Lmin to Lmax- We denote the number density 
of galaxies in the range of luminosity considered as n, in units of Mpc~^. Thus, the number density of 
galaxies per unit luminosity is ^{L) = n$(L). For plotting purposes, we use the luminosity function per 
unit magnitude, normalized to the mean density of galaxies: 

$(M) = 0.41n(10)nL$(L), (4) 

with M = —2.5\ogiQ{L/Lo), where Lq is the luminosity of an object with absolute magnitude zero. 

We use the maximum likelihood approach outlined by Efstathiou, Ellis, & Peterson (1988) and Sandage, 
Tammann, & Yahil (1979) to calculate ^{L). This method is based on the conditional probability of 
observing a galaxy with a luminosity between Lj and Lj + dLj given its redshift Zj: 

^' ' p{zj)dzj J^--^^[^^ dL^L) ' 

Here, L^i^^zj) and Lyna,x{zj) are the minimum and maximum luminosities observable at redshift zj, given 
the faint and bright flux limits of the survey. As described below, we evaluate Equation (5) using two models 
for $(i): a Schechter function and a non-parametric model. Note that the probability in Equation (5) is 
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completely independent of the density field of galaxies; this fact makes the method insensitive to the effects 
of large-scale structure, at least to the extent that the luminosity function is universal. In addition, this 
insensitivity makes it necessary to calculate the normalization in a separate step, described below. 

The procedure is to minimize the quantity 

C = -2Y,Wo\og\p{Lj\zj)], (6) 

i 

where the weight Wc is applied to properly account for galaxies which are eliminated due to fiber collisions. 
That is, we assign double weight to a galaxy if it has eliminated a neighbor from the redshift sample due to 
a collision. In principle, we should use the weight wd ft, to also weight those galaxies that do not have fibers 
due to being in dense regions, but the changes in our results would be extremely small. Note that the weights 
we assign only affect the mean of the maximum likelihood estimator if they correlate with the luminosity; 
that is, the weights Wc only make a difference to the extent that the typical luminosities of galaxies involved 
in fiber collisions are different from those of other galaxies. 

To perform the minimization, we need to choose a model for Here we use both a Schechter (1976) 

function and a non-parametric model. The Schechter function is 

$(L)dL = (^^^ exp dL. (7) 

Expressed per unit magnitude, it is 

$(M)dM = 0.41n(10)(/.aO-°-^(^-*^-)("+^) exp LlQ-O-^^^-^-)] dM. (8) 

In this case, the parameter space of and a is searched to minimize C. To calculate the error bars in this 
case, WC! perform 200 Monte Carlo calculations of the following form. We take all the galaxy positions as 
observed in the survey, and for each assign a random luminosity drawn from the fitted luminosity function 
and the fiux limits. Then we fit for the best Schechter function in each Monte Carlo sample and use the 
variation among samples to calculate the errors and covariance of M* and a. 

We have checked the difference between the Monte Carlo error estimate and a jackknife error estimate 
(Lupton 1993). For the jackknife method, we divide our sample into eighteen regions on the sky of approx- 
imately equal area. Then we perform the luminosity function analysis eighteen times, each time leaving a 
different section out. The estimated statistical variance of a parameter x is 



WB.r{x) = ^^Y.^x-xf, (9) 

i=l 

where TV = 18 in this case and x is the mean value of the parameter measured in the samples. The errors 
derived in this way are within 30% of those found using the Monte Carlo method described in the previous 
paragraph. 

The second model for $(i) we use is the non- parametric model, which is a piece- wise constant interpola- 
tion in logarithmically spaced steps in luminosity, as described by Efstathiou, Ellis, & Peterson (1988). For 
this model, the likelihood given by Equation (5) can be maximized quickly by using an iterative equation. In 
this case, because there are so many parameters to constrain in the covariance matrix, we follow Efstathiou, 
Ellis, & Peterson (1988) and calculate the errors based on the second derivatives of the likelihood function 
around the solution. 
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Given the fit to citlicr model, wc need to calculate the normalization. We do so here using the minimum 
variance estimator of Davis & Huchra (1982): 

" JdV^{z)w{z) ' ^ ^ 

where the integral is over the volume covered by the survey between the minimum and maximum redshifts 
used for our estimate. The weight for each galaxy is 

and the selection function is 

/■imax(z) 

0(0)= / dL^L)ft, (12) 

where ft is the galaxy sampling rate as described above. The integral of the correlation function is: 



/>oo 

J3= / drr^^ir). 
Jo 



(13) 



Clearly, because n appears in the weight w{z), it must be determined iteratively, which we do using the 
simple estimator n — {1/V) '^wd <i>{zj) as an initial guess. We iterate until the change is less than 10% of 
the estimated error. The error in the normalization is given by 



J dV(t){z)w{z) 



-,1/2 

(14) 



We have also checked the error estimates in the normalization using the jackknife method described above 
and find the results to be consistent with those of Equation (14). In addition, there are hints in the jackknife 
results that the errors in n may be correlated with the errors in M* and a. However, eighteen jackknife 

samples still provide a rather noisy estimate of the six parameters in a 3 x 3 covariance matrix. We leave 
the question of these correlated errors for a future analysis of a larger data set. 

The luminosities of galaxies are known to correlate with other quantities, such as their surface brightness, 
color, and morphology. We investigate these relationships below and will need a way to calculate the joint 
relationship between luminosity and other quantities. To do so, we use a two-dimensional generalization 
of the method of Efstathiou, Ellis, & Peterson (1988), along the lines of Sodrc & Lahav (1993), de Jong 
& Lacey (2000), and O'Neil & Bothun (1999) (all of whom considered the specific case of galaxy surface 
brightnesses). Instead of bins in luminosity, we take bins in the two-dimensional plane of luminosity and the 
other quantity in question. The derivation of the method of Efstathiou, Ellis, & Peterson (1988) and its 
implementation are virtually unchanged from the one- dimensional case. 

The redshift and flux limits we apply in each band are tabulated in Table 1. The high redshift limit 
is imposed to minimize noise in the determination of the normalization. The low redshift limit is applied 
to avoid difficulties in the commissioning data associated with the accidental deblending of nearby, large 
galaxies. While the flux limit in r* is set by the spectroscopic flux limit, we calculate the luminosity function 
in other bands by defining flux limits in these bands sufficiently bright that almost all of the galaxies brighter 
than the limit were actually targeted. Imposing the requirement that < 2% of galaxies brighter than the 
flux limit in the given band are fainter than r* = 17.6 yields the flux limits of wfjj^ — 18.4, g^:^^ = 17.65, 
**im — 16.9, and z*^^ = 16.5. That is, for u* and g* we make sure that the flux limit is bright enough to 
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include all but the very bluest galaxies at r* = 17.6, and for i* and z* we make sure it is bright enough 
to include all but the very reddest galaxies. This approach guarantees that the full range of galaxy colors 
is represented at each apparent magnitude in each sample. Neglecting the dependence of the i^-corrections 
on galaxy color, this would be the equivalent of considering nearly the full range of galaxy colors at each 
rcdshift. Note that for u* and g*, the flux limits arc shallow enough that we need to reduce the maximum 
redshift we consider. To use a larger fraction of the data to obtain the same result, we could calculate the 
bivariate luminosity function in r* and in the other band of interest, but for the moment we will stick to the 
simple and conservative approach. 



Using the galaxy sample described in Section 2, we fit the models for the luminosity function described in 
Section 3. The results for the r* band are shown in Figure 5 for the case of the flm = 0.3, Ha — 0.7 cosmology. 
The parameters of the fit to the Schechter function for all three models are given in Table 2, along with 
their errors and covariances. The covariance is expressed as the correlation coefficient r between the errors. 
As previous authors have found, the errors in and a have a strong correlation. Judging from the non- 
parametric results, the Schechter function is a good fit to the data. In fact, from the likelihoods calculated 
in the fits, we found that the Schechter function is a better fit than the non-parametric fit. This result is 
almost certainly due to the fact that the non-parametric method assumes that the luminosity function is flat 
within each bin; for example, if we double the number of luminosity bins for the non-parametric flt, then its 
likelihood is higher than that of the Schechter function fit. Future work will investigate this issue further 
using non-parametric methods like those of Koranyi & Strauss (1997) and Springel & White (1998), which 
adopt more general interpolation schemes than the standard piece-wise constant interpolation of Efstathiou, 
Ellis, & Peterson (1988). 

We show the results of our non-parametric fits and Schechter function fits for each of the four other 
bands in Figure 6 and Table 2. The Schechter function fit is good in every band; its slope a for low luminosity 
galaxies is somewhat steeper in u*, but in general it is a remarkably weak function of bandpass. The non- 
parametric fits are well-constrained down to absolute magnitudes of —16 to —17.5, depending on bandpass, 
about four magnitudes less luminous than M* in the central {g*, r* , i*) bands. 

If we extrapolate the fitted Schechter function to low luminosities, the luminosity density can be calcu- 
lated: 



Similarly, one can add up the luminosity in the non-parametric fit over the absolute magnitude range that we 
have measured it. We list in Table 2 the luminosity density j found for the extrapolated Schechter function 
in each band (in units of absolute magnitude per hr^ Mpc'^) and the fraction of this luminosity /„,, that is 
accounted for in the non-parametric fit. The errors given in j are calculated from the errors in a, M.^,, and 
, accounting properly for the covariance between the errors of a and but assuming the errors in ^* are 
independent; these errors arc consistent with the errors estimated based on the jackknife samples described 
in Section 3. From these results, we see that if the fitted slope a continues to describe the luminosity function 



4. Results for the Luminosity Function 



4.1. Luminosity Function in Five Bands 




= (l)M{a + 2). 



(15) 
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at low luminosities, the galaxies represented in our sample account for most of the luminosity density of the 
universe. We will discuss the caveats associated with surface brightness selection effects later in Section 4.3. 

Table 2 also lists the luminosity density expressed in solar units. To determine these quantities, we 
assume Mvq = 4.83, {B - V)q = 0.65, {U - B)q = 0.13, and {R - /)© = 0.34 (Binney & Merrifield 1998), 
and use the stellar color transformations of Fukugita et al. (1996) to obtain 

M^.Q = 6.39 ; Mg.Q = 5.07 ; M,.© = 4.62 ; M^.© = 4.52 ; M,.© = 4.48. (16) 

These numbers agree well with those predicted from theoretical models of G2V stars (Lcnz et al. 1998: Fan 
1999). There still may be errors at the level of 5% in the photometric calibration for the data analyzed here; 
furthermore, the color transformations are somewhat uncertain, perhaps also by 5%. For this reason we add 
in quadrature a 7% systematic uncertainty to the errors listed for the luminosity density in solar units. The 
luminosity densities determined with this analysis are consistent with those determined by Yasuda et al. 
(in preparation) based on the bright number counts in a similar SDSS data set along with the parameters 
M* and a of the Schechter function measured here, providing a check on our method of determining the 
normalization. 



4.2. Redshift Distribution of SDSS Galaxies 

Here we compare the observed redshift distribution to the expected redshift distribution of a homo- 
geneous galaxy sample based on our luminosity function and the flux limits. We calculate the expected 
number of galaxies using a "typical" galaxy color of g* — r* = 0.65 in order to calculate the ii'-corrcction. 
One can refine this analysis using the joint color-luminosity relation we calculate below in Section 4.4, and 
it is important to do so when calculating large-scale structure statistics. However, for current purposes 
picking a typical color will suflice. To calculate the observed redshift distribution we weight each galaxy 
according to its "collision weight," Wc, described in Section 2.6. Figure 7 compares the expected counts to 
the observed counts. Clearly large-scale structure is noticeable in an area of sky this small, but overall the 
redshift distribution appears consistent with that predicted by the luminosity function. Zehavi et al. (in 
preparation) arc; measuring the large-scale structure in this sample and find that it is generally consistent 
with what has been found in previous surveys. We should emphasize again that the likelihood method used 
here to calculate the shape of the luminosity function is insensitive to the effects of large-scale structure, to 
the extent that the luminosity function does not depend on environment. 

To reassure ourselves that the method of targeting galaxies does not introduce redshift-dependent se- 
lection, apart from the ciTcct of the flux limits, we divide the set of galaxies into two subsamples above and 
below cz = 30,000 km s^^, and calculate the luminosity function separately for these two samples. Figure 
8 compares the results. The two luminosity functions agree in the absolute magnitude range in which they 
overlap. More precise comparisons of the normalization of luminosity functions at high and low redshift will 
be possible once more area of the sky is available for analysis. 

4.3. Surface Brightness Selection Effects and the Correlation of Luminosity with Surface 

Brightness 

Our assertion that we are measuring the galax;y luminosity function implicitly assumes that our sample 
is truly flux-limited, and that surface brightness selection effects do not influence the luminosity function at a 
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significant level. As described in Section 2.4, the commissioning observations incorporate an explicit surface 
brightness cut at Hi/2.r* = 23.5. In the apparent magnitude range 14.5 < r* < 17.6 used here, we expect 
the SDSS imaging to detect essentially all galaxies above this surface brightness threshold, so there should 
be no additional selection effects associated with failure to detect low surface brightness galaxies in the first 
place. The half-light radius of a galaxy with apparent magnitude m and half-light surface brightness /ii/2 
is ri/2 = 6.0" X ioO-2(/'i/2-23.5)-K0.2(i7.6-m)_ jv^gg^^. ^.j^g magnitude limit r* = 17.6 and the surface brightness 
limit /ii/2 = 23.5, this size is well below the ~ 50" scale used to define the background for sky subtraction. 
At r* = 14.5, the half-light radius is 25", so sky subtraction and deblending failures could have some 
impact on detection, but only a tiny fraction of the galaxy targets woiild be influenced, many of which would 
be too close in redshift to be included in the sample used here for the luminosity function (Yasuda et at, in 
preparation, discuss this issue in more detail). At the depth of the SDSS imaging a galaxy with r* = 17.6 and 
A*i/2,r* = 23.5 has mean S/N ^ 50 within ri/2, so detections have high statistical significance. In principle, 
one can also miss high surface brightness galaxies by mis-classifying them as stars, but Figure 3 shows that 
the distribution of the quantity mpsp — TOmodei used for star-galaxy separation has fallen off well before 
reaching the classification threshold. The number of high surface brightness galaxies missed is therefore a 
tiny fraction of the full sample; if wc choose a higher cut-off for mpsF — "^modoi which excludes 1% of the 
canonical galaxy sample, the effect on the luminosity and surface-brightness distributions is negligible. 

The /ii/2,r* = 23.5 surface brightness limit itself only excludes about 0.5% of galaxies within our flux 
limits. However, the well-known relationship between surface brightness and luminosity (Lilly et al. 1998; 
Driver 1999; de Jong & Lacey 2000), confirmed below, suggests that the excluded galaxies will preferentially 
be of low luminosity. Since these underluminous galaxies are observable only within a rather small, nearby 
volume, they are given disproportionate statistical weight when the luminosity function is calculated, and 
thus even a small number of missing objects may significantly bias our estimate of the luminosity function 
at low luminosities. This effect is one reason why the full SDSS sample will have a fainter surface brightness 
limit than 23.5, as mentioned in Section 2.4 above. 

A simple way to investigate this effect is to raise the surface brightness threshold to Ati/2,r* = 22.5. This 
restriction excludes an additional 1.5% of the galaxies within our flux limits. Figure 9 shows the luminosity 
function of the full and restricted samples in the top panel, and the ratio of the luminosity functions in 
the bottom panel. Significant numbers of galaxies are missed at Mr* > —19 in the restricted sample. In a 
Schcchtcr function fit to the samples, the slope a. is shallower by about 0.08 in the restricted sample, and 
the influence on a is approximately the same in the other bands. For galaxies in the absolute magnitude 
range —23.5 < Mr* < —16, the restricted sample has about 5% lower luminosity density than does the full 
sample, as measured from the non-parametric fits. If one extrapolates the fitted Schechter functions to all 
luminosities, the difference is still only 5%. However, if the low-luminosity slope a steepens considerably 
below M > —16, then the contribution of low surface brightness galaxies to the luminosity could be larger. 
Because the objects with /Ui/2,r* > 23.5 have a surface number density that is about 0.5% of that of galaxies 
in our sample, these effects represent a firm upper limit on the influence of our surface brightness threshold 
on the luminosity function in the range of absolute magnitudes considered here. 

The joint distribution of surface brightness and luminosity is a valuable diagnostic for understanding 
surface brightness selection effects and an important quantitative characterization of the galaxy population 
in its own right (see, e.g. Driver & Cross 2000 and references therein). To define this distribution, we 
take /ii/2.r* (corrected for cosmological surface brightness dimming) as the measure of surface brightness, 
and use the two-dimensional version of the non-parametric method of Efstathiou, Ellis, & Peterson (1988) 
to calculate the joint distribution of /Xi/2,r* and Mr*, taking into account the apparent magnitude and 
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surface brightness limits of the sample. The results of this fit are shown in Figure 10. Figures 11 and 12 
show slices parallel to the axes of this two-dimensional plane: the luminosity function for several ranges of 
intrinsic surface brightness (Figure 11) and the surface brightness distribution for several ranges of absolute 
magnitude (Figure 12). Note that we have not corrected the values of /Ui/2,r* for the effects of seeing, so the 
surface brightness of smaller galaxies is systematically underestimated. 

Figures 10 — 12 demonstrate a number of interesting points. Galaxies with M^, < — 19 show a clear peak 
in their surface brightness distribution, with typical /ii/2.r* ~ 20. This result is in rough agreement with 
the previous results of Freeman (1970) and Courteau (1996), although those authors measured the central, 
rather than the half-light, surface brightness, and concerned themselves only with disk galaxies. In addition, 
galaxy surface brightness is clearly correlated with luminosity, in agreement with the results of, e.g., dc Jong 
& Lacey (2000). For lower luminosities, the typical surface brightness is lower, and the distribution of surface 
brightness is broader (Phillipps & Disney 1986). For lower surface brightness, the galaxy luminosity function 
becomes steeper. The intrinsic surface brightness distribution of galaxies with M^.* < — 19 drops off well 
before our threshold of /ii/2.r* = 23.5, providing further evidence that our luminosity function is unaffected 
by surface brightness selection in this magnitude range. For galaxies less luminous than Mr- « —17.5, the 
surface brightness distribution does decline before hitting our threshold, but not in as definitive a manner as 
it does for more luminous galaxies; thus, our galaxy density in this regime might be underestimated, and we 
cannot rule out an upturn in the luminosity function at these magnitudes. However, the luminosity density 
for a Schechter luminosity function with a > —2 is dominated by galaxies with M « M*. Unless there is a 
radical change in the shape of the luminosity function at Mr» > —17.5, we can safely say that low surface 
brightness galaxies make a relatively small contribution to the luminosity density of the universe. All of 
these conclusions are in broad agreement with most recent studies (Impey & Bothun 1997). 

There is great potential for further studies of galaxy surface brightness and its correlation with other 
properties using the SDSS. As a result of improvements in the photometric data reduction software, the 
main survey (now underway) will adopt a fainter surface brightness threshold, H\/2,r* = 24.5. This change 
will improve the measurement of the luminosity function and the bivariate luminosity-surface brightness 
distribution at lower luminosities. With the aid of simulated data, it should be possible to correct estimated 
/ii/2 values for seeing effects. With a larger data sample, it will be possible to analyze joint distributions of 
surface brightness with a larger set of other galaxy parameters, such as color and morphology, using all five 
bandpasses. Finally, the SDSS imaging will provide an excellent data set for finding candidate galaxies far 
below the Hi/2,r-' = 24.5 threshold, even if their redshifts are not obtained by the SDSS spectroscopic survey. 
Follow-up efforts in optical and 21 cm and correlation with HI radio surveys such as HIPASS {e.g., Banks 
et al. 1999) will greatly improve the characterization of the low surface brightness galaxy population. 

4.4. Dependence of Luminosity on Color 

It is well-known that galaxies of different intrinsic colors have quite different luminosity functions. We 
explore this issue here in a simple way, by calculating the joint distribution of luminosity and the intrinsic 
g* — r* color. We determine the intrinsic g* — r* color for a galaxy at a particular redshift by using the 
observed g* — r* color to interpolate between the predicted colors at that redshift for various galaxy types in 
the results of Fukugita, Shimasaku, & Ichikawa (1995). There are only five templates between which we are 
interpolating, which may affect the accuracy of our color estimates (typical color corrections to rest frame 
are about 0.25 magnitudes, so our accuracy is considerably better than that). 
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The resulting distribution, $(Mr»,(.g* — r*)o) is shown in Figure 13. This figure shows the "E/SO" 
ridge at (<?* — r*)o ~ 0.75, which is well-known in the context of galaxy clusters, as well as the dependence 
of color on luminosity. As in the case of surface brightness, we show the luminosity function for several 
ranges of color in Figure 14, and the color distribution as a function of luminosity in Figure 15. The ridge 
at {g* — r*)o ^ 0.75 is apparent at high luminosities in Figure 15, as is the relative flatness of the color 
distribution at low luminosity. 

A better examination of the dependence of galaxy color on luminosity will be possible once the photo- 
metric calibration of the data is complete (although the calibrations used here are thought to be accurate to 
0.05 magnitudes). Once the final calibration is available, another improvement will be to use the broad-band 
imaging (and possibly the spectra) to determine the fC-corrections from the data itself. Methods developed 
for photometric redshifts, which can use the broad-band colors to reconstruct the typical spectral energy 
distributions of galaxies {e.g., Csabai et al. 2000), can be brought to bear on this problem and will resolve 
the difficulties mentioned above associated with the small number of templates used for interpolation. 

4.5. Dependence of Luminosity on Morphology 

Finally, let us consider the dependence of luminosity on the morphology of galaxies. Although there 
are numerous ways to quantify galaxy morphology, we choose here for simplicity to use the "concentration 
index," defined to be c = rgo/r^o, the ratio of the radii containing 90% and 50% of the Petrosian flux. We 
will actually study the inverse concentration index 1/c here, because it has the property of being between 
zero and one. The concentration index is high (c « 3.3, or 1/c « 0.3) for pure de Vaucouleurs profile galaxies 
and low (c « 2.3, or 1/c « 0.43) for pure exponential profiles. It has been shown to correlate well with 
visual morphological classification for bright galaxies (Shimasaku et al., in preparation; Strateva et al, in 
preparation). 

Figure 16 shows the joint distribution of luminosity and inverse concentration index. As one might 
expect, the high luminosity galaxies tend to have high concentration index, and thus to have profiles which 
are closer to de Vaucouleurs profiles. As before, we slice this two-dimensional plane to get the luminosity 
function for several ranges of the inverse concentration index in Figure 17, and the distribution of the inverse 
concentration index for several ranges of luminosity in Figure 18. 

Bear in mind that the concentration index is not the only, and perhaps not the best, measure of 
morphology. For example, some of the inverse concentration indices measured may be artificially affected by 
the effects of seeing. The inverse concentration index for a realistic PSF measured from the survey is around 
1/c ~ 0.5, so exponential disks and de Vaucouleurs profiles that are marginally resolved will tend to have 
their inverse concentration indices increased towards this value. Furthermore, this measure of morphology is 
based only on the radial profile in r* , and thus does not use the information on color gradients, shape, and 
galaxy substructure that is available in the imaging. Developing measures of morphology that use the full 
information from the galaxy images and that are independent of the effects of seeing is one of the ongoing 
projects within the SDSS collaboration. 



As mentioned in Section 2.3, rso and rgo axe defined with respect to the Petrosian flux; the inverse concentration parameter 
for de Vaucouleurs profiles when rso and rgo are measured with respect to the total fiux is about 1/c 0.18. 
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5. Comparison of Petrosian Magnitudes to Isophotal Magnitudes 

As discussed in Section 2.4, we liave selected spectroscopic targets and calculated galaxy luminosities 
using Petrosian magnitudes, so that the fraction of a galaxy's light that is measured is independent of the 
amplitude of its surface brightness profile and independent of cosmological redshift dimming or Galactic 
extinction. In order to compare our luminosity function results to those of previous surveys, most of which 
selected galaxies according to some variant of an isophotal magnitude, we must understand how the difference 
between Petrosian and isophotal magnitudes affects the determination of the luminosity function. 

We can calculate an isophotal magnitude for each galaxy in our sample, by integrating the azimuthally- 
averaged radial profile of each galaxy, which is available as a standard output of photo (Lupton et al, 
in preparation), out to a chosen isophote. These magnitudes do not correspond exactly to the common 
definition of isophotal magnitudes, which typically are not azimuthally-averaged. However, as we show in 
the next section by comparing the SDSS results to those of other surveys, the two definitions are probably 
nearly equivalent, at least for luminosity statistics. 

Before selecting galaxies according to our newly defined isophotal magnitudes and calculating a lumi- 
nosity function, let us first directly examine the relationship between the SDSS Petrosian magnitudes and 
isophotal magnitudes. If mjso — Wpetro were the same for all galaxies, the luminosity scale would differ 
in the zero-point, but there would be; no other change in the shape of the luminosity function. However, 
rriiso — TOpetro depends on both surface brightness (because low surface brightness objects hide a higher frac- 
tion of their light below any given isophotal level) and redshift (because of cosmological surface brightness 
dimming); thus the resulting luminosity function is strongly affected. To show the redshift dependence. 
Figure 19 shows rriiso ~ 'Tipotro for a volume-limited sample of galaxies to z = 0.2 (roughly Mr* < —21.7), 
for several choices of the isophotal limit. In each panel we show a linear regression fit, showing clearly that 
the isophotal flux is not a constant fraction of the Petrosian flux as a function of redshift. Figure 20 shows 
miso — TOpetro as a function of absolute magnitude. For luminous galaxies, rriiso — Wpetro is large because 
most such galaxies are near the edges of the survey (because more volume is there) and thus suffer consider- 
able surface brightness dimming. For underluminous galaxies, Wiso — Wpetro is again large, now because (as 
dictated by the correlation of luminosity and surface brightness found in Section 4.3) many of these objects 
arc; low surface brightness. The redshift and luminosity dependence persists even at isophotal limits so faint 
that the typical redshift zero isophotal magnitude is brighter than the same galaxy's Petrosian magnitude. 
For isophotal limits more typical of other surveys, such as /Xr*,iso — 23, which is comparable to that of the 
LCRS, or at /Ur*,iso = 24, which under the assumption that {bj — r*) w 1 is comparable to the 2dFGRS, the 
trend with redshift is quite strong. In Section 6, we will perform a much more careful comparison of our 
results to those of the LCRS and the 2dFGRS. 

Figure 21 compares the luminosity functions derived using the isophotal and the Petrosian samples in 
the r* band. Clearly there is a considerable difference between the estimates, which is not attributable 
to a simple constant offset between isophotal and Petrosian magnitudes. Instead, the isophotal magnitude 
estimate differs from the Petrosian magnitude estimate at the low and high luminosity ends. Figure 21 also 
lists /lum, the fraction of the integrated luminosity density in the Petrosian sample that is recovered in each 
isophotal sample. Prom these results, we infer that surveys with shallower isophotal limits could be missing 
considerable amounts of luminosity density. We will compare our results more directly to other surveys (and 
in the appropriate bands) in Section 6 below. 

There are two reasons that isophotal samples can underestimate the luminosity density when the isopho- 
tal limits are too bright. First, isophotal magnitudes with bright limits measure less flux for each galaxy, as 
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revealed by Figure 20. Second, because the fraction of the total flux measured by an isophotal magnitude 
decreases with redshift from cosmological surface brightness dimming, using the standard formula for the 
distance modulus overestimates the effective volume of each galaxy (Dalcanton 1998). The dominant effect 
in Figure 21 is the first, that less light is measured. To show this, we recalculate luminosity functions using 
a simple 1/Vmax estimator (e.g., Binggcli, Sandagc & Tammann 1988). Wc first use the Pctrosian magni- 
tudes both to determine the absolute magnitude of the galaxy and to apply the flux limits (and thus also 
to determine Vmax)- Then we use the isophotal magnitudes to do both, finding the expected reduction in 
luminosity density (about 40% in the case of /Ur'.iim = 23). We then calculate the luminosity function again, 
determining Mr- from the isophotal magnitude but applying the flux limits to the Pctrosian magnitude. In 
this case, the reduction of the luminosity density is nearly as great, around 35%. Finally, we determine M^. 
from the Petrosian magnitude, but apply the flux limits to the isophotal magnitude, thus isolating the effect 
of an incorrect VJnax determination. In the case that /ir'.iim = 23, this procedure reduced the luminosity 
density only by 10%. These results show that it is the missing light beyond the isophotes that dominates the 
reduction in the luminosity density, rather than an overestimate of the maximum volume. Presumably, this 
result is related to that of Dalcanton (1998), that the systematic effects associated with isophotal magni- 
tudes have only a moderate effect on the estimated V/Vmax distribution. The error in V^nax estimated with 
the standard distance modulus may be large in some individual cases, but most galaxies enter the sample 
at a redshift close to the largest one at which it is possible for them to, because that is where most of the 
volume is. The additional dimming out to the galaxy's limiting redshift is therefore small in most cases. 



6. Comparison to Other Surveys 

As wc have seen above, the choice of isophotal magnitudes can greatly affect the results for the galaxy 
luminosity function, confirming the results of Dalcanton (1998). Thus, in order to show that the results of 
the SDSS are consistent with those of other surveys, we must reanalyze the SDSS data in the same ways that 
previous surveys analyzed their data. The SDSS provides the ability to do so (for surveys of comparable; 
depth) in a way that no other large-scale survey has before. Using the five-band color information and 
the color conversions for galaxies provided by Fukugita, Shimasaku, & Ichikawa (1995), one can convert 
the SDSS magnitudes into the appropriate band for almost any desired optical survey. We convert to the 
LCRS R band and to &j, and compare to the LCRS and 2dFGRS luminosity fimctions. finding significantly 
more luminosity density in both cases. Using the measured photometric properties, one can reconstruct an 
isophotal magnitude similar to that of any given survey (albeit with azimuthally averaged light profiles). 
When we do so for the LCRS and for 2dFGRS, we find our results are consistent with theirs. These results 
demonstrate that the LCRS and 2dFGRS are missing significant fractions of the light in the universe because 
their measurements of galaxy magnitudes are based on fairly bright isophotes. 



6.1. Comparison to the LCRS 

The LCRS photometry is expressed in its own variant of the R band, a hybrid of Gunn r and the 

Kron-Cousins R band, which wc will denote i?GKC following Shectman et al. (1996). This band can be 
related to SDSS magnitudes as follows (Fukugita, Shimasaku, & Ichikawa 1995; Shectman et al. 1996): 



i^GKC =r* - 0.05 - 0.089(5* - ^*)- 



(17) 
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We use this transformation to calculate i?GKC-band Petrosian magnitudes for all SDSS galaxies, corrected 
for Galactic extinction. Using the method described in Section 3 for defining complete samples in u*, g*, i* 
and z*, we define a complete sample in i?GKC by setting a fltix limit of mum = 17.40. For consistency with 
Lin et al. (1996), who calculate the iicKC-band luminosity function in the LCRS, we use an Einstein-de 
Sitter universe in this case. The results are listed in Table 3. The top panel of Figure 22 shows the R- 
band luminosity function determined from the SDSS in this way and compares our result to that of Lin et 
al. (1996). The two luminosity functions diverge at low and high luminosities for the reasons discussed in 
Section 5: isophotal magnitudes underestimate the luminosities of low luminosity galaxies because of their 
low intrinsic surface brightnesses, and of high luminosity galaxies because of cosmological surface brightness 
dimming. The ratio of the i?-band luminosity density found by Lin et al. (1996) to that found here is 
/lum = 0.48. That is, there is (at least) twice as much i?GKC-band light in the universe than found by 
the LCRS. Table 3 lists these results, along with the conversions of the luminosity densities to solar units 
(assuming Mj^^-^^q = 4.52, using the results of Binney & Merrifield 1998 and the relation i?GKC = -Rkc+0.1 
from Shectman et al. 1996; these choices are consistent with the results quoted by Lin et al. 1996). 

Let us now analyze the SDSS images using the same method as used in the LCRS. First, we use 
isophotal magnitudes limited at /iijcKcHm = 23; this number approximately corresponds to the statement 
of Shectman et al. (1996) that the LCRS isophotal magnitudes are limited at 15% of the sky brightness. 
Second, Shectman et al. (1996) exclude galaxies of low "central surface brightness," as measured by the 
magnitude within a fixed angular aperture about the size of a fiber. Similarly, we apply a "central surface 
brightness" cut based on the "fiber magnitude," an aperture magnitude with a 3" diameter. Following 
Shectman et al. (1996), we define the fiber magnitude cut to be: 

miiGKcfiber < 18.85 - 0.5(17.7 - tur^^^), (18) 

Finally, for consistency we calculate the i^-corrections using the formula K{z) = 2.51ogio(l + z), as Lin et 
al. (1996) did (although this makes rather little difference in our results). The bottom panel of Figure 22 
shows the comparison of the resulting luminosity hmction to that of Lin et al. (1996). They are largely 
consistent, except that the SDSS sample has a slightly brighter M*; this difi'erence may persist because our 
isophotal limit is not precisely what LCRS used. Nevertheless, the general consistency shows that the extra 
fiux detected in the SDSS in the top panel is most likely real, and that the missing light in the LCRS is 
due to having shallow isophotes and excluding galaxies with low central surface brightnesses. To understand 
how much is due to the central surface brightness cut of Equation (18), we list results in Table 3 for the 
case in which we do not apply this cut; this test indicates that the central magnitude cut is responsible for 
a substantial portion of the underestimate in the luminosity density. 



6.2. Comparison to the 2dFGRS 

The 2dFGRS sample is based on APM plates (Maddox, Efstathiou, & Sutherland 1990) and uses the 
bj band. Here we relate the SDSS photometry to bj using the relation 

6j =5* + 0.14 + 0.088(5* -r*), (19) 

determined using the results of Fukugita, Shimasaku, & Ichikawa (1995) and the galaxy color relation 
bj = B — 0.35(5 — V) from Metcalfe, Fong, & Shanks (1995). We follow the same procedure as for the 
LCRS in the previous section to create a sample limited at mbj^n-m = 17.8. Again, for consistency with 
Folkes et al. (1999), we assume an Einstein-de Sitter universe. While our results are much closer to those of 
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2dFGRS than they were to the LCRS, we still measure about 1.4 times the luminosity density that Folkes et 
al. (1999) do, as we show in Figure 23. This difference is significant, given that the stated errors in the bj 
luminosity density in the 2dFGRS are about 5%, and the errors for the SDSS are about 15%. Table 3 lists 
these results, now using M^.q = 5.30 to convert to solar units (determined using the stellar color relation 
bj = B - 0.28(S - V) from Blair & Gilmore 1982 and the results of Binney & Merrifield 1998). 

Let us analyze the SDSS images using the methods used by 2dFGRS. APM magnitudes are isophotally 
limited at ^i,. = 25 mag per square arcsecond, and then corrected to total magnitudes by assuming that the 
galaxy profile is a circularly symmetric Gaussian (Maddox et al. 1990). After determining these "APM- like" 
magnitudes for SDSS galaxies, and applying a flux limit of mbj^nm = 17.8, we again calculate the luminosity 
function. This flux limit is considerably shallower than that of the 2dFGRS. Thus, the mean rcdshift is 
lower, and to accurately match the isophotal limits would require a slightly brighter isophotal limit; we have 
instead decided to err on the side of including a little more flux. The result is shown in the bottom panel 
of Figure 23. The SDSS luminosity function determined in this manner is nearly identical to the 2dFGRS 
luminosity function. (Although the Schcchtcr parameters arc rather different, it is clear from the flgurc that 
the values are nearly degenerate). Given our estimated errors and the fact that our bj sample is considerably 
shallower than that of the 2dFGRS, this agreement of the two results is somewhat surprising. This result 
indicates that the 40% difference in luminosity density between the SDSS and 2dFGRS is due to the fact 
that the APM isophotal limits exclude a significant amount of light in the outer regions of galaxies, even 
when corrected. For comparison, we also give results in Table 3 for the case of bj = 25 isophotal magnitudes, 
with no correction to "total" magnitudes based on the Gaussian model. These results show that while the 
APM corrections work in the right direction, the Gaussian model is an insufficient description of observed 
galaxy profiles for extrapolation to "total" magnitudes. 



7. Conclusions 

We have presented the luminosity function and luminosity density of galaxies in five optical bandpasses 
(m*, g*, r* , i* and z*) using 11,275 galaxies in SDSS commissioning data. The Schechter function appears 
to be a good fit in all cases, with a low-luminosity slope that is remarkably similar in all bandpasses. On 
the other hand, the slope a is sensitive to galaxy surface brightness, color, and morphology. The redshift 
distribution of the galaxies is generally consistent with our derived luminosity function, although large-scale 
structure is still evident in our sample. 

We have demonstrated the correlation between r* luminosity and galaxy surface brightness, color, and 
morphology. In particular, we have shown that luminous galaxies tend to be higher surface brightness, 
redder, and more concentrated than less luminous galaxies, in accordance with previous results. We find 
that galaxies obey a strong magnitude-surface brightness relation. There is a high characteristic half-light 
surface brightness for luminous galaxies; at the same time, low-luminosity galaxies show a broad surface 
brightness distribution, with many low surface brightness objects. As described above, refinements of the 
photometric analysis beyond the standard software used for the survey will improve our understanding of 
these correlations. For example, a better understanding of the effects of seeing, inclination, and internal dust 
extinction will improve our measures of surface brightness and morphology. Improvements in our estimates 
of /sT-corrections, which will be possible using both the spectroscopy and the five-band photometry of the 
SDSS, will improve our understanding of the distribution of galaxy colors and allow us to estimate the stellar 
mass function of galaxies. Work in progress on the spectral classification of SDSS galaxies {e.g., Castander 
et al. 2001) is showing the same trends of galaxy type with luminosity found here. 
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We have shown that a considerable fraction of the luminosity density in the universe has been missed 
by previous samples that were based on shallower imaging. We measure 2 times the i?GKC-band luminosity 
density found in the LCRS and 1.4 times the 6j-band luminosity density found in the 2dFGRS. Work in 
progress, which directly compares the measured magnitudes of galaxies in these surveys to SDSS magnitudes 
in the regions where the surveys overlap, will test whether these difFcrcnccs arc indeed reasonable. While 
large-scale structure could still have some effect on results from these surveys, our internal tests show that 
the differences of results are almost certainly due to our adoption of Petrosian magnitudes in preference to 
isophotal magnitudes. In particular, we can reproduce the LCRS and 2dFGRS results if we mimic their 
photometry and selection effects. The high luminosity density found by the SDSS implies that the cosmic 
density of stellar matter is higher than previously thought. If this is the case, it may relieve some of the 
strain on galaxy formation models, which typically predict a higher density of stellar matter than previous 
observational estimates (Katz, Weinberg, & Hernquist 1996; Pearce et al. 1999; Granato et al. 2000). 

Our estimates of the luminosity density, given in Tables 2 and 3, are almost certainly still underestimates. 
First, for the case of dc Vaucoulcurs profiles, we know that up to 18% of the galaxy light may be missed by the 
Petrosian magnitudes as we define them. Including this light would increase somewhat the total luminosity 
density found here (though only by a few percent, because de Vaucouleurs profile galaxies constitute only 
about 20% of all galaxies). Second, we have shown that while the explicit surface brightness limit /Ui/2,iim = 
23.5 of the commissioning data exclude only 0.5% of galaxies, these galaxies very likely represent a larger 
fraction of the luminosity density than that; for example, only 1.5% of galaxies are in the range 22.5 < 
Mi/2,r* < 23.5, but such galaxies are responsible for ~ 5% of the luminosity in the universe. The deeper 
surface-brightness limit of the main SDSS survey will thus provide a more complete estimate of the luminosity 
density in the universe. Third, our estimates assume a constant low luminosity slope a. An upturn in 
the luminosity function at absolute magnitudes less luminous than M = —16, suggested by some studies 
(Loveday 1997; Phillipps et al. 1998), would mean there was yet more luminosity not included in this 
analysis. A final possibility, related to the previous one, is that a significant fraction of stars exist outside of 
galactic environments, having perhaps been tidally stripped from their galaxy of birth. As more SDSS data 
is obtained and we perform more sophisticated analyses of the various selection effects in the survey, we will 
be able to address many of these problems. 

Tables of the non-parametric fits to the luminosity function in each band and the joint relationships 

between luminosity and surface brightness, color, and concentration, are available on the World Wide Wcb^^. 
We expect that these results and the parametric fits in Table 2, which represent a straightforward analysis 
of only a small fraction (one percent) of the expected SDSS sample, will already provide interesting new 
constraints on theories of galaxy formation. In addition, this characterization of the luminosities, surface 
brightnesses, colors, and morphologies of local galaxies provides a solid baseline for interpreting the evolution 
of galaxies observed at higher redshift. Future work with the SDSS will allow us to explore the correlations 
between all of these properties and to describe fully the population of galaxies in the local universe. 
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Table 1. Limits Applied to Sample in Each SDSS Band 



Band 


Flux Limits 


Redshift Limits (km s ^) 


Number of Galaxies 


u* 


14.50 < m < 18.40 


5,000 <cz< 30,000 


1,679 


9* 


14.50 < TO < 17.65 


5, 000 <cz< 50, 000 


4,684 


r* 


14.50 < m < 17.60 


5,000 <cz< 60,000 


11,275 


i* 


14.50 < m < 16.90 


5,000 <cz < 60,000 


7,441 


z* 


14.50 < m < 16.50 


5,000 <cz< 60,000 


6,090 



Table 2. Parameters of Fits to Luminosity Function in SDSS 







Band 


4>. 

(10-2 Mpc-3) 


M* - 5 logio h 




a 


rM,,a 


j (mags) 


j {h W^Lq) 


/np 


0.3 


0.7 


u* 


4.00 ±0.90 


-18.34 ±0.08 


-1 


.35 ±0.09 


0.80 


-15.21 ±0.26 


4.35 ± 1.08 


0.85 






9* 


2.06 ±0.23 


-20.04 ± 0.04 


-1 


.26 ±0.05 


0.79 


-16.05 ±0.13 


2.81 ±0.38 


0.94 






r* 


1.46 ±0.12 


-20.83 ±0.03 


-1 


.20 ± 0.03 


0.78 


-16.41 ±0.09 


2.58 ±0.28 


1.00 






i* 


1.28 ±0.11 


-21.26 ±0.04 


-1 


.25 ±0.04 


0.77 


-16.74 ±0.10 


3.19 ±0.37 


0.98 






z* 


1.27 ±0.11 


-21.55 ±0.04 


-1 


.24 ±0.05 


0.74 


-17.02 ±0.11 


3.99 ±0.48 


0.97 


0.3 


0.0 


u* 


4.48 ± 1.06 


-18.27 ±0.07 


-1 


.32 ±0.09 


0.78 


-15.20 ±0.26 


4.34 ± 1.07 


0.84 






9* 


2.48 ± 0.33 


-19.97 ±0.05 


-1 


.23 ±0.05 


0.83 


-16.16 ±0.15 


3.09 ±0.47 


0.98 






r* 


1.71 ±0.15 


-20.74 ±0.03 


-1 


.20 ±0.03 


0.74 


-16.49 ±0.09 


2.78 ±0.31 


0.96 






i* 


1.47 ±0.13 


-21.17 ±0.04 


-1 


.23 ±0.04 


0.77 


-16.78 ±0.10 


3.31 ±0.39 


0.96 






z* 


1.45 ±0.14 


-21.46 ±0.04 


-1 


.20 ±0.05 


0.80 


-17.02 ±0.11 


3.99 ±0.48 


1.00 


1.0 


0.0 


u* 


4.68 ±1.13 


-18.24 ±0.07 


-1 


.31 ± 0.09 


0.82 


-15.20 ±0.26 


4.34 ± 1.10 


0.84 






9* 


2.42 ±0.30 


-19.92 ±0.04 


-1 


.22 ±0.05 


0.80 


-16.07±0.13 


2.87 ±0.40 


0.97 






r* 


1.87 ±0.18 


-20.67 ±0.03 


-1 


.15 ±0.03 


0.84 


-16.47 ±0.10 


2.74 ±0.31 


1.00 






i* 


1.61 ± 0.15 


-21.11 ±0.04 


-1 


.20 ±0.05 


0.78 


-16.78 ±0.11 


3.32 ± 0.40 


0.98 






z* 


1.62 ±0.16 


-21.39 ±0.04 


-1 


.19 ±0.05 


0.74 


-17.08 ±0.11 


4.19 ±0.51 


0.96 



Note. — Luminosity function in various bands assuming the three cosmological models considered here, with Hq = 100 
km s~^. The Schechter function is characterized by the three parameters ^* (in units of 10-^ Mpc-^), M*, and 
a, whose values and errors in each case are listed here. Also listed is the covariance between the errors of M* and a, 
determined from the Monte Carlo error determination, j is the luminosity density as determined by integrating the 
Schechter function over all luminosities (again in units of h"^ Mpc-"^). We list j both in absolute magnitudes and in 
solar luminosities (determined as described in the text). The error bars on the value of j in solar luminosities include 
a 7% contribution due to calibration and color transformation uncertainties, /np is the fraction of the luminosity in 
the extrapolated Schechter function which is accounted for in the non-parametric fit; that is, this is the fraction of the 
inferred light which comes from galaxies with luminosities which are observed in the sample. Naturally, it is lowest when 
the slope a at low luminosities is highest. 
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Table 3. Comparison of SDSS Luminosity Function to 2dFGRS and LCRS 



Band 


Sample 


4>* 


M* — 5 logio 


a 


j (mags) 








(10-2 MpC-3) 










^G;K-C 


SDSS Petrosian 


1.92 ±0.23 


-20.80 ±0.03 


-1.17±0.03 


-16.64 ±0.12 


2.91 ±0.38 




SDSS IsophotaP 


2.13 ±0.26 


-20.37 ±0.04 


-0.94 ±0.05 


-16.15 ±0.12 


1.86 ±0.24 




SDSS Isophotal'' 


2.02 ± 0.24 


-20.30 ±0.04 


-0.74 ±0.05 


-15.95 ±0.12 


1.54 ±0.20 




LCRS 


1.90 ±0.1 


-20.29 ±0.02 


-0.70 ±0.03 


-15.87±0.1 


1.4±0.1 


bj 


SDSS Petrosian 


2.69 ±0.34 


-19.70 ±0.04 


-1.22 ±0.05 


-15.97 ±0.14 


3.21 ±0.46 




SDSS Isophotal'= 


2.17 ±0.28 


-19.55 ±0.05 


-1.12 ±0.06 


-15.47 ±0.14 


2.04 ±0.30 




SDSS Isophotal'i 


2.32 ±0.30 


-19.62 ±0.04 


-1.15 ±0.05 


-15.65 ±0.13 


2.39 ± 0.34 




2dFGRS 


1.69 ±0.17 


-19.73 ±0.06 


-1.28 ±0.05 


-15.56 ±0.05 


2.19 ±0.12 



Note. — Same as Table 2, for the LCRS i?G;K-c-band and the bj band. We include comparisons to the LCRS 
(Lin et al. 1996) and the 2dFGRS (Folkes et al. 1999). 

"^Without the limits on mfiber in Equation (18). 

^With the limits on mfiber in Equation (18). 

'^Without APM-like corrections to the isophotal magnitudes. 

^With APM-like corrections to the isophotal magnitudes. 
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Fig. 1. — Behavior of Petrosian magnitudes as a function of galaxy size, for pure, face-on exponential profiles 
(top panel) and circularly symmetric de Vaucouleurs profiles (bottom panel), under marginal seeing conditions 
in the SDSS. For each profile, the fraction of total light measured by the Petrosian flux is plotted against the 
(true) angular half-light radius of the galaxy. In the absence of seeing, Petrosian flux measures a constant 
fraction of the total light as a function of redshift, unlike isophotal fluxes. However, small differences do 
appear as a function of the angular size due to seeing. As the galaxy size becomes comparable to the size of 
the seeing disk, the fraction of light measured by Petrosian quantities approaches that fraction for a PSF, 
which is about 95%. In the case of exponential disks, this reduces the flux, because in the absence of seeing 
nearly 100% of the light is measured. In the case of a de Vaucouleurs proflle, this increases the flux, because 
in the absence of seeing about 80% of the light is measured (assuming that the de Vaucouleurs profile extends 
to infinity). As shown in Figure 2, a significant number of galaxies have sizes close to that of the seeing disk, 
so this dependence on galaxy size for small galaxies (which is present in isophotal magnitudes as well) will 
limit the accuracy of our estimates of the total luminosity density, though only to a few percent. 
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Fig. 2. — Distribution of half-light radii rso of SDSS galaxies with r* < 17.6. Top panel shows the distri- 
bution; bottom panel shows the fraetional cumulative distribution. Roughly half of the galaxies are in the 
regime (rso < 2.5") in which seeing significantly affects the determination of the Petrosian radius and flux. 
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Fig. 3. — Distribution of the quantity wpsp — Wmodei (determined in r*), which is used to separate stars from 
galaxies for the purposes of target selection in the SDSS. The cut-off of 0.242 used in the commissioning data 
for this purpose is shown as the vertical line. The solid histogram is the distribution of topsf — Wmodei for 
targeted objects which were spectroscopically identified as galaxies; the dotted histogram is the distribution 
for targeted objects which were spectroscopically identified as stars. About 1% of objects targeted as galaxies 
turned out to be stars. The distribution falls off sharply at small mpsF — Jrimodeh implying that <C 1% of 
galaxies brighter than r* = 17.6 are misclassified as stars. 




Fig. 4. — ii'-corrections in the r* band for galaxies with intrinsic g* — r* colors of 0.2, 0.4, 0.6, and 0.8, 
as labeled. These curves are obtained by linear interpolation between the morphological types used by 
Fukugita, Shimasaku, & Ichikawa (1995). Similar results are available for the other four bands. 
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Fig. 5. — Top panel shows the luminosity function of SDSS galaxies in the r* band, based on a sample of 
11, 275 galaxies in the commissioning data. For these results, we assumed Sim = 0.3 and Oa = 0.7. Results 
for other cosmologies are given in Table 2. The Schechter fit {thick line) and the non-parametric fit {thin 
line with error bars) agree well. The correlation coefficient between a and is r = 0.78, showing that the 
errors in these parameters are strongly correlated. Bottom panel shows the number of galaxies in each bin 
used to determine the luminosity function; even though the least luminous galaxies have the highest space 
density, there are very few such galaxies in a flux-limitcd sample. The fact that the low-luminosity end of 
the correlation function is based on a relatively small fraction of the galaxies in the sample is the reason that 
the errors in the normalization of the luminosity function are not comparable to the Poisson errors expected 
given the number of galaxies; that is, a small number of the galaxies are given disproportionate statistical 
weight, making the errors greater than Poisson. 
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Fig. 6. — Same as Figure 5, for the u* , g* , i* , and z* bands. The flux hmits and number of galaxies in each 
sample arc listed in Table 2. As in the r* luminosity function in Figure 5, the Schechter fits {thick line) and 
the non-parametric fits {thin line with error bars) agree well. 
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Fig. 7. Rcdshift distribution of SDSS galaxies (histogram) compared to the average distribution expected 
given the r* luminosity function and the flux limits (using if-corrections for a "typical" galaxy color of 
g* — r* = 0.65). The redshift distribution is reasonable, although large-scale structure is obvious. 
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Fig. 8. — Comparison of the luminosity function in r* determined using samples above and below cz = 30, 000 

km s^^. In the region in which the two luminosity functions overlap they are in good agreement, indicating 
that the redshift dependence of the selection function is purely due to the flux limits of the survey. More 
precise tests of this nature can be made when larger areas of the sky (which are less aff^ected by large-scale 
structure) become available. 
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Fig. 9. — Luminosity function of SDSS galaxies in the r* band for the full sample, limited at a surface- 
brightness of /ii/2.r* < 23.5, and for a restricted sample with /ii/2,r* < 22.5. The higher surface brightness 
threshold excludes about 1.5% of the galaxies in the full sample. The top panel shows both luminosity 
functions, showing that differences between them appear at Mr- > —19. The bottom panel shows the ratio 
of the two luminosity functions. 
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Fig. 10. — The logarithmic grcyscalc represents ^{Mr* , fii/2.r*) in for SDSS galaxies, in units of galaxies 
per Mpc'^ per unit magnitude per unit surface brightness, calculated using the method of Efstathiou, 
Ellis, & Peterson (1988). The translation of the greyscale to these units is given in the bar at the top of 
the figure. The data sample has an explicit half-light surface brightness cut at /Ui/2,r* = 23.5. The strong 
correlation between surface brightness and luminosity is apparent. It is clear why many of the objects with 
A*i/2,r* > 22.5 are low luminosity, as shown in Figure 9. 
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Fig. 11. — Luminosity function of SDSS galaxies in the r* band, in bins of half-liglit surface brightness, in 
units of galaxies per Mpc"^ per unit magnitude. Each curve shows the luminosity function for a range of 
surface-brightnesses. As labeled, the curves are colored blue, cyan, green, yellow, and red from high to low 
surface brightness. The steep faint-end slope for low surface-brightness galaxies is apparent, as well as the 
fact that the high surface-brightness galaxies tend to be luminous. 
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Fig. 12. — Half-light surface-brightness function of SDSS galaxies in the r* band, in bins of galaxy luminosity, 
in units of galaxies per Mpc'^ per unit surface brightness. Each curve shows the surface brightness 
function for a range of Mr' . As labeled, the curves are colored blue, cyan, green, yellow, and red from low 
to high luminosity. The luminous galaxies tend to have a characteristic half-light surface brightness, while 
low luminosity galaxies have a broad distribution that extends to low surface-brightness. 
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Fig. 13. — Same as Figure 10, for the joint luminosity-color distribution ^(Mr*,^* — r*) in r* for SDSS 

galaxies, now in units of galaxies per Mpc^ per unit magnitude per unit color. The intrinsic g* —r* colors 
are inferred using the observed g* — r* color, the measured redshift, and the results of Fukugita, Shimasaku, 
& Ichikawa (1995). The "E/SO" ridge at — r* « 0.75 is apparent, as is the strong correlation between 
luminosity and color. 



- 40 - 




Fig. 14. — Same as Figure 11, only now showing the luminosity function in bins of intrinsic g* — r* color, 
in units of galaxies per Mpc'^ per unit magnitude. As labeled, the curves are colored blue, cyan, green, 
yellow, and red from bluest to redder. The intrinsic color is determined using the observed color, the 
measured redshift, and the results of Fukugita, Shimasaku, & Ichikawa (1995). Errors in the calibration 
used for this data probably broaden this distribution slightly, by about 0.05 magnitudes. 
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Fig. 15. — Same as Figure 12, only now showing the distribution of intrinsic g* ~ r* color for several ranges 
of r* luminosity, in units of galaxies per h^^ Mpc"^ per unit color. As labeled, the curves are colored blue, 
cyan, green, yellow, and red from low to high luminosity. As in Figure 13, the "E/SO" ridge at g* —r* « 0.75 
is apparent for high luminosities, while for low luminosities, the g* — r* distribution is nearly flat. 
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Fig. 16. — Same as Figure 10, for the joint distribution of luminosity and inverse concentration index 
$(Mr*, 1/c) in r* for SDSS galaxies, in units of galaxies per Mpc"^ per unit magnitude per unit inverse 
concentration index (ICI). The concentration index is defined to be c = rgo/rso- Again, the galaxies with a 
low inverse concentration index {i.e., those galaxies with profiles closer to de Vaucouleurs) tend to be brighter. 
Note, however, the abundance of high inverse concentration index galaxies (1/c > 0.45). Since exponential 
disks should have inverse concentration indices around 0.43, these galaxies either have much broader profiles 
than exponential, or are being affected by seeing and noise in the measurement. Nevertheless, the trend of 
luminosity with morphology is clear. 
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Fig. 17. — Same as Figure 11, only now showing the luminosity function in bins of inverse concentration 
index 1/c = ?'5o/'"90) in units of galaxies per Mpc^ per unit magnitude. As labeled, the curves are 
colored blue, cyan, green, yellow, red, and black from least to most concentrated. 
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Fig. 18. — Same as Figure 12, only now showing the distribution of inverse concentration index 1/c — r^o/rgo 
for several ranges of r* luminosity, in units of galaxies per h^^ Mpc'^ per unit inverse concentration index. 
As labeled, the curves are colored blue, cyan, green, yellow, and red from low to high luminosity. 
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Fig. 19. — Comparison of isophotal and Petrosian magnitudes as a function of redshift in the r* band for 
several choices of the isophotal limit. We only show galaxies luminous enough to be included in the sample 
out to z = 0.2 {Mr* < —21.7). We have quantified the dependence on redshift with a linear regression along 
z. 
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Fig. 20. — Comparison of isophotal and Pctrosian magnitudes as a function of absolute magnitude in the 
r* band for several choices of the isophotal limit. We show all galaxies with r* < 17.6. The most luminous 
galaxies have high miso — tripetio because they tend to be further away, and thus suffer relatively more surface 
brightness dimming. The least luminous galaxies have high rriiso — ^ripetro because they tend to be low surface 
brightness. 
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Fig. 21. — Fitted r* luminosity function for Petrosian magnitudes (upper curve in each panel, same as Figure 
5) compared to the luminosity fiuiction using isophotal magnitudes (lower curve in each panel). Each panel 
corresponds to a different choice of limiting isophotal surface-brightness fir* Aim, as labeled. The Schcchter 
function parameters for each choice of isophotal limit are listed in each panel, as well as the fraction /lum of 
the total Petrosian luminosity that is detected in each isophotal sample. 



-48- 



0.1 k 
^ 10-2 

J 10-3 

10-5 



0.1 = 



^ 10-2 

o 

I 

<^ 10-4 
10-5 



G;K-C 



-band Sample 



SDSS R^_^^ ^ 
"(Petrosi'an) 

= 1.92 ± 0.23 (xlO-2 
M^= -20.80 ± 0.03 
a = -1.17 ± 0.03 




- SDSS R 



LCRS Sample (Lm, et al. 1996^ 
fium = 0-48 

0^= 1.90 ± 0.10 (xlO-2) 
^ -20.29 ± 0.02 
= -0.70 ± 0.05 



G;K-C 



band Sample 
^Isopholal with ytij- =23) 
= 2.02 ± 0.24 TxlO-s) 



_M^= -20.30 ± 0.04 



= a 



-0.74 ± 0.05 




LCRS Sample (Lm, et al. 1996) 



1.90 ± 0.10 (xlO-2) 
-20.29 ± 0.02 
-p. 70 ± 0.05 



24 



-22 



20 



18 



16 



M 



R, 



G:K-C 



Fig. 22. — i?-band luminosity function in the SDSS compared to that in the LCRS calculated by Lin et 
al. (1996). Top panel shows the luminosity function calculated from the SDSS using Pctrosian magnitudes. 
The ratio of the i?-band luminosity density in the LCRS to that in the SDSS is /lum = 0.48. The bottom panel 
shows the results of analyzing the SDSS data in the same manner as the LCRS, using isophotal magnitudes 
limited at /UH,iim = 23 and applying a central magnitude cut similar to that of the LCRS, as described in the 
text. This sample is considerably more similar to the results of the LCRS itself. The remaining discrepancies 
may be due to our use of a slightly deeper isophotal magnitude than that of the LCRS, the axisymmetric 
nature of our isophotal magnitudes, or the slightly different definition of our central magnitude cut. 
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Fig. 23. — Same as Figure 22, this time comparing the SDSS bj luminosity function to that of the 2dFGRS as 
calculated by Folkes et al. (1999). The top panel shows that the SDSS sample based on Petrosian magnitudes 
has about 1.4 times the luminosity density of the 2dFGRS luminosity function. The bottom panel shows 
that if we use isophotal magnitudes with the same isophotal limit that the 2dFGRS uses, and "correct" these 
magnitudes in the same manner as they do, we obtain nearly the same results. These results indicate that 
the 2dFGRS misses a significant amount of luminosity density due to its relatively shallow isophotal limits, 
notwithstanding the efforts to correct the isophotal magnitudes assuming a universal Gaussian profile for 
galaxies. 



